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Abstract. As some of the first known objects to exist in the Universe, Lyman 
alpha emitting galaxies (LAEs) naturally draw a lot of interest. First discovered 
over a decade ago, they have allowed us to probe the early Universe, as their 
strong emission line compensates for their faint continuum light. While initially 
thought to be indicative of the first galaxies forming in the Universe, recent 
studies have shown them to be increasingly complex, as some fraction appear 
evolved, and many LAEs appear to be dusty, which one would not expect from 
primordial galaxies. Presently, much interest resides in discovering not only the 
highest redshift galaxies to constrain theories of reionization, but also pushing 
closer to home, as previous ground-based studies have only found LAEs at z > 3 
due to observational limitations. In this review talk I will cover everything from 
the first theoretical predictions of LAEs, to their future prospects for study, 
including the HETDEX survey here in Texas. 



1. Introduction 

As new and more powerful instruments and telescopes are developed, we now 
probe deeper into the Universe than ever before. While once redshifts of a 
few tenths were difficult to observe, we can now routinely probe to z = 4-7, 
corresponding to a time when the Universe was only ~ 1 Gyr old. The goal is 
simple in principle. Whe n we look at t he galaxies around us, we see a well defined 
morphological sequence ( Hubbldll926l ). extending out to z > 1. We would like to 



understand how the first galaxies formed, and through what evolutionary means 
they evolved into the galaxies we see today. 

To study the early Universe, we need techniques to accurately determine the 
distances to galaxies. Locally, various distance ladder rungs are used, but the 
most powerful of these, SNe la, taps out at z ~ 1. At high redshift, a common 
method is to use the expected spectral fe atures of the galaxies we h ope to find, 
such as with the Lyman break technique ( Steidel fc Hamilton 19931 ). At 2: > 2, 



the Lyman limit, an intrinsic spectral break at 912 A, shifts into the observed 
optical. Though the intergalactic medium (IGM) is largely ionized, enough rem- 
nant neutral gas is in a typical line-of-sight to efficiently absorb photons bluer 
than the Lyman limit, increasing the amplitude of this break. Additionally, 
from 912-1216 A, a forest of absorption lines appear due to absorption of inter- 
vening Hi along the line of sight, known as the Lya forest. At higher redshifts, 
the number of absorbers increases, and the opacity in the Lya forest increases 
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dramatically, such that the spectral break begins at 1216 A, and is called the 
Lya break. 

Using this spectral feature, one can discover large samples of Lyman break 
galaxies (LBGs) by observing with three broadband filters. For example, at 
z 3, the Lyman break is in the n'-band, thus one would observe with u' , 
g' and r' filters. As the break passes through the u'-band, the u' — g' color 
will increase. By examining how the colors of stellar population model tracks 
evolve with redshift, one can select a. u' — g' color which will select galaxies at 
2.5 < z < 3.5. However, this single color can also include red galaxies at lower 
redshifts, thus one typically also demands a blue rest-frame ultraviolet (UV) 
color, g' — r' in our example. With this method, one can efficiently select large 
numbers of high-redshift star-forming galaxie^. 

However, the LEG selection method relies on robust continuum detections 
in multiple images, as well as a very deep dropout filter. Thus, typic al LEG stud- 
ies p lace a brightness limit on their sample, typically R < 25.5 ( Steidel et al.l 
19961 ). leaving a large population of faint galaxies unexplored. Additionally, 
though galaxies selected with this method should lie in a redshift slice of Az 1, 
one requires followup spectroscopy to confirm the selection criteria and deter- 
mine the redshift selection function. This can be costly, as even "bright" LEGs 
(m-AB ~ 23) can require many hours on 8-lOm class telescopes to confirm their 
redshifts via absorption lines. 

If the goal is to study the earliest galaxies, then one can rely on some addi- 
tional spectral features. In particular, they sho uld be heavily star-forming , and 
thus bright in numerous nebular emission lines. Partridge fc Peebles! ( 19671 ) pre- 
dicted that one could select the first galaxies on the basis of their Lya emission, 
which at a rest wavelength of 1216 A shifts into the optical from 2 < z < 6. With 
a typical initial mass function (IMF), many stars will form with sufficient mass 
(and thus temperature) to emit large quantities of hydrogen-ionizing photons 
(A < 912 A). These photons wih ionize Hi out to approximately the Stromgren 
radius. At this point, there is a balance between ionization and recombination, 
and when a proton an d electron recomb in e, a Lya photon will be emit ted ~ 67% 
of the time (Case E, IOsterbock|[l989l 1. iPartridge fc Peebles! (|l967l ) predicted 
that > 5% of the total luminosity of these galaxies could emit in this single line. 

Spurred by these predictions, a number of studies sea rched for these Ly- 
man alpha emitting gal axies (LAEs), but to no avail (e.g., Koo fc Kron!ll980! : 
Meier fc Terlevichl 198l! ). A number of causes were suggested to explain the 
lack of LAEs, such as t he possibility that these galaxies could be dusty (e.g. 



Meier &: Terlevich 198ll ). In a uniform interstellar medium (ISM), Lya photons 



will resonantly scatter off Hi, resulting in longer path lengths, and thus they 
stand a much greater chance of being absorbed by dust. Or, some reasoned that 
perhaps these objects simply didn't exist, or were far too faint to be observed. 

T hankfully, for tunes turned in the late 19 90's, as two groups searching for 
LAEs, [C owie &: Hul (|l998l ) and iRhoads et all (|200Cll ). discovered large samples. 
Why were they able to do this while previ ous groups w e re no t? The answer 
is a matter of technologica l advancemeii t. Cowie &: Hu ( 19981 ) used the new 
Keck 10m telescope, while iRhoads et al.l (j2000! ) used the new MOSAIC large 
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format detector on the Kitt Peak National Observatory 4m telescope. The 
deeper imaging afforded by Keck or the wider field-of-view provided by MOSAIC 
both resulted in a larger volume being probed, and thus gave these studies a 
greater chance of discovery. 

These, and subsequent studies, used the method of narrowband selection 
to identify these galaxies on the basis of their bright Lya lines. To do this, 
one observes with a narrowband filter, where the redshift of the LAEs one 
find s is defined by the c entral wavelength and width of the filter. For exam- 
ple, iRhoads et al. ( 2000l ) used a set of filters around 6500 A, thus the LAEs 
they found were at z ~ 4.5. One also observes the field with a broadband filter 
encompassing the narrowband filter. Examining sources in both images, one se- 
lects objects with a flux excess in the narrowband image as possible emission line 
candidates. However, low-redshift interlopers are possible; for example, with a 
narrowband filter at 6500 A, one could also select [On], [Oiii] or H/3 emitters. 
To screen against these, one typically enforces a dropout criterion similar to the 
Lyman break method, requiring an additional broadband image blueward of the 

line. 

Using this method of narrowband selection. ICowie &: Hul l 19981 ) and lRhoads et al 



tOOd ) both found that one could discover large samples of LAEs fairly easily, 
at least at redshifts where Lya is observed in the optical. One of the first inter- 
esting discoveries about LAEs was that their Lya equivalent widths (EWs) were 
strong (the EW is a measure of the line strength relative to the continuum). 
In many cas es, the EWs exceede d predictions from normal stellar population 
models (e.g.. ICharlot fc Falll[l993l ). Figure [J shows the distribution of Lya EW 
versus age for different star- formation histories (SFHs). Unless one is observing 
a LAE in the first few Myr, the Lya EW should be < 200 A. However, this is 
frequently not the case (e.g., iKudritzki et"all l2000l : iMalhotra fc Rhoadd [200I ) . 
Thus, something unexpected is occurring in many LAEs. 



2. Uncovering the Physical Properties of LAEs 

There are a number of physical possibilities which could result in these high EWs. 
If LAEs are truly the first galaxies, or their immediate descendants, we might 
expect them to have very low, or perhaps zero metallicities. The lack of metal 
lines for cooling during star formation can thus result in a n initial mass fu nction 
(IMF) which preferentially forms more massive stars than Salpeter ( 19551 ). or a 



so-called top-heavy IMF. Heavier and more metal-poor stars have hotter stellar 
photospheres, which result in the production of more ionizing photons, and 
thus more Lya photons. To examine the possibility of this scenario will require 
detailed stellar population modeling or direct measurements of the metallicity 
from rest-frame optical emission lines. 

If these galaxies host super-massive black holes which are actively accreting 
material, then much of the light we see could originate from the accretion disk 
around the black hole, including any Lya emission. The Lya EW would thus not 
be subject to such strict limits as in the case of star formation. If many LAEs 
hosted active galactic nuclei (AGN), this could solve the EW problem. However, 
at high redshift the only means we currently have to diagnose the presence of 
AGN is via their X-ray emission. Typically LAE studies examine X-ray data. 
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Figure 1. Model variation of Lya EW with age for three different star 
formation histories: constant, exponentially decaying, and a simple stellar 
population, or burst (from Finkelstein et al. 2008). All of these models have 
Z = 0.02 Zq, and include the effect of IGM absorption at z ~ 4.5, where we 
assume that half of the line flux (i.e. the blue wing) will be subject to IGM 
attenuation. 



and they find that the AGN contamination is low, from 1-5% at 2 < z < 5 (e.^ 
Malhotra et"al1l2003l : IWang et al.ll2004l : iGawiser et alll2006l : iNilsson et aDl200£ 



Thus, AGN are likely not a major contaminant, although low-lumi nosity AGN 
may b e lurking beneath current X-ray sensitivities (see, e.g., iFinkelstein et all 
2009bl ). 



A different scenario for the high EWs observed in LAEs involves a dusty 
ISM. While dust was originally thought of as the reason why LAEs were not 
seen in the 1980's, this relied on the assumption of a homogeneous ISM. If, 
however, the ISM was clumpy, where neutral gas and dust were mixed together 
in clouds with a relatively empty inte r-cloud inediu r n, it is possible to a ctually 
enhance the Lya we would measure (|Neufeldl[T99ll : iHansen Ohll2006l 1. This 
can happen because Lya photons are resonantly scattered by Hi, while contin- 
uum (and other emission lines) are not. Continuum photons traveling through 
such an ISM would penetrate the surface of these clouds. Some of these photons 
will travel straight through and escape the galaxy, some will be absorbed by 
dust, and some will be scattered by the dust into a different direction. How- 
ever, Lya photons will resonantly scatter with the Hi on the surface of these 
clumps, essentially bouncing off. Thus, Lya photons stand a significant chance 
of evading interaction with dust in such an ISM. This scenario has the net ef- 
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feet of increasing the ratio of the Lya-to-continuum escape fractions, and thus 
increasing the observed EW. 

2.1. Spectral Energy Distribution Fitting 

The typical method of uncovering the physical properties of distant galaxies 
is through spectral energy distribution (SED) fitting, whereby one compares 
observations of a galaxy (in m ost cases broadband photometry) to a s uite of 
stellar population models (e.g.. iPapovich et al.ll200l1 : IShaplev et al.|[200ll ). One 



varies the age, dust extinction, metallicity, SFH, etc. of the models, and finds 
the one that best represents the observations, typically by minimizing the 
between the models and the observations. 

The first SED analysis of LAEs hoped to answer the question of whether or 
not they truly were the primitive first galaxies. However, many of these studies 
relied on ground-based data, thus they had to stack the fluxes of their sample 
together t o obtain the necessary signal-to-noise (S/N) to robustly compare them 
to models. iGawiser et al.l ( 20061 ) studied 18 spectroscopically confirmed LAEs at 
z = 3.1, finding that on avera ge they had ages o f 90 M yr, masses of 5 x 10^ Mq 



z = 6.L, imdmg that on avera ge they haa ages o i 90M yr, masses or b x iO Mq, 
and zero dust extinction. In iFinkelstein et al.l ( 20071 ), we stacked nearly 100 



LAEs at z ~ 4.5 into three bins by EW strength, and found age and mass ranges 
from 1-40 Myr, and 10''"^ Mq, respectively, with the highest EW sources having 
the youngest ages and lowest masses. 

In the f i rst st udies to individually analyze LAEs, iLai et af] ( 20071 ) and 



Pirzkal et al.l ( 20071 ) found strikingly differe nt results, althou gh the differences 



can be explained by their selection biases. iLai et all ( 20071 ) studied LAEs at 



z ~ 5.7 in the Great Observatories Origins Deep Survey (GOODS) North field, 
restricting their analysis to three LAEs which were detected in the infrared 
(IR) with the Spitzer Space Telescope. They found that these objects had ages 
from 5-700 Myr, masses from 5 x 10^^^*^ Mq, and significant dust extinction 
from Ay = 0.4-1.7 mag. However, the restriction to Spitzer detections likely 
biases towards more massive LAEs, though t he existence o f LAEs this mas- 
sive is intriguing. In a complementary study, iPirzkal et al.l ( 20071 ) studied 10 



LAEs individually at 4 < 2 < 5.5 discovered in the Hubble Ultra Deep Field 
(HUDF) via grism spectroscopy with the Advanced Camera for Surveys (ACS) 
on board the Hubble Space Telescope {HST). This method of discovery in the 
extremely small HUDF will likely bias thi s study towards int rinsically fainter 
and lower mass objects. Correspondingly, iPirzkal et al. ( 20071 ) find ages of 1- 



20 Myr, M = 3x 10^-2 x 10^ Mq and extinction of Ay = 0-0.6 mag. Thus, the 
faintest LAEs are among the lowest mass and youngest galaxies in the Universe. 

2.2. Investigating the ISM Geometry 

However, these early studies did not use an important piece of information we 
have on these objects, namely that they exhibit Lya emission. Lya can have a 
significant effect on the observed SED , increasing the broa dband magnitude by 
up to 0.5 mag for strong emitters (e.g.. iPapovich et al.ll200ll ). thus Lya emission 
should be included in the models. Additionally, since the majority of LAEs 
were discovered via narrowband photometry, the narrowb and data should be 



included in the fit. However, most studies use the popular iBruzual Sz Chariot 
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( 20031 . hereafter BC03) stellar population synthesis models, which do not include 
nebular emission lines. 



NB665/673 ACS/F435W (B) F606W (V) F775W (i') F850LP (z') IRAC/3.6|ijm 4.5>im S.B^ni S-Ourn 




Figure 2. 10" cutouts of seven of the z ~ 4.5 LAEs from (from Finkelstein 
et al. 2009a). They are extremely bright in the narrowband and undetected 
in the B-band. 11/14 LAEs were detected in at least one Spitzer band. 



This was the motiva tion behind our work in Finkelstein et al. ( 20081 ) and 



Finkelstein et al.l (l2009el h We obtained narrowband imaging from the Cerro 



Tololo InterAmerican Observatory 4m telescope, using the MOSAIC II wide- 
field imager in three filters: Ha (Ac = 6563 A), Ha+80 (Ac = 6650 A) and [Sii] 
(Ac = 6725 A). These filters can thus select LAEs at 4.37 < z < 4.57. We chose to 
observe the Extended Chandra Deep Field South (ECDF-S), which contains the 
GOODS-South field, consisting of extremely deep HST and Spitzer data. We 
found 14 LAEs in the GOODS region, restricting our analysis only to objects in 
the deep HST fields so that we could analyze these LAEs individually. Figure [2] 
shows image cutouts of seven of our LAEs from this study. 

To take advantage of the measured Lya properties of these objects (from 
their narrowband excesses), we added Lya emission to the models, basing the 
line strength off the number of ionizing photons at a given age, metallicity and 
SFH. Adding i n this emission also allowed us to examine the dust enhancement 
scenario (e.g., iNeufeldl Il99ll ) discussed earlier, as we were able to attenuate 



Lya by a different amount than the continuum when dust extinction was added 
to the models. Our treatment was slightly more complex than simple clumpy 
spheres, as we allowed the amount that dust affected Lya for a given model to 
be a free parameter, which we called "g", where q is defined by: 

fLya = fhja * exp{-qTc) (1) 
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where Tc is the optical depth of dust to the continuum around Lya. A value of 
q = thus symbolizes the Neufeld scenario, where Lya escapes any dust atten- 
uation, and q = 10 symbolizes a homogeneous ISM, where resonant scattering 
ensures that Lya will be attenuated much more than the continuum. 

When all of this was included in our models, we found the interesting result 
that LAEs are not a homogeneous population of galaxies. Out of our sample 
of 14 galaxies, we found that 12 appeared to be young, with best-fit ages of 
3-15 Myr and masses of 0.16-6 X W^Mq. However, the remaining two LAEs 
were much older, both with ages of ~ 500 Myr, and masses of 0.4-5 x 10^^ Mq. 
Perhaps more surprising is that all of the galaxies in our sample appear to be 
dusty, with Ay ~ 0.1-1.1 mag. However, the best-fit dumpiness parameters 
(q) show that the ISM may in fact be enhancing the observed Lya EWs, as 
10/14 objects had best-fit g-values of < 1, implying the EWs we observe have 
been enhanced over the EWs intrinsic to the star- forming regions. We thus 
concluded that ISM enhancement of the Lya EW is a plausible explanation for 
the observed plethora of high-EW sources. Additionally, the ubiquity of dust 
implies that LAEs are not the primitive first galaxies as previously thought. 
Finally, while the LAE population appears to be predominantly young, there is 
a non-negligible fraction of older, more evolved sources. The relation between 
these rarer, old LAEs and the younger, more common LAEs is an important 
topic for future research. 



3. What the Near Future Holds 



The previous section has illustrated that at the very least, LAEs are a com- 
plicated category of objects. The cause of the large observed EWs among 
narrowband-selected LAEs is still not clear, though we have shown that EW 
enhancement via a clumpy ISM is consistent with some observations, and thus 
is a plausible scenario. However, it is clear that on average, LAEs are young 
(< 100 Myr) a nd low mass (< 10 ^ Mf?)). Additionally, studies of LAEs from 
3 < z < 6 re.g.. lGawiser et alJl2006l : iPirzkal et aIll2007l : lFinkelstein et aDl2009el ) 
show that their physical properties do not appear to change with redshift, includ- 
ing their physical s ize (Malhotra et al. 20 10, in prep). This is contra ry to LBGs 
which get larger (Ferguson et al. 20041 ) and more massive (e.g., Stark et al 



20091 : iFinkelstein et alJl2009dl ~ with decreasing redshift. This suggests a sce- 



nario where, at every redshift, LAEs are the build ing blocks of more evolved 
galaxies at lower redshifts (e.g.. ICawiser et al.l 120071 1. 

In order to test this theory and truly understand the physical make-up 
of LAEs, we need to directly measure their properties. For example, dust 
extinction is one of the least constrained properties deduced from SED fit- 
ting. Fortunately, the UV energy that the dust absorbs is emitted in the 
rest-frame far-infrared (FIR). Unfortunately, at high-redshift this emission is 
shifted into the sub-millimeter-millimeter regime. In the future, the Atacama 
Large Millimeter Array (ALMA) should be able to dete ct the dust emission from 
^ 50% of LAEs in a relatively sh ort amount of time ([Finkelstein et al.ll2009d : 
Daval. Hirashita. Sz Ferrarall2009l ). but this is years away. 

Currently, near-infrared (NIR) spectroscopy can resolve some of these is- 
sues. For example, the Balmer decrement provides a direct measurement of the 



8 



Finkelstein 



dust extinction. For this experiment, one measures H a and H/3 emiss ion. In 
the absence of dust, the ratio of Ha/H/3 should be 2.86 ( Qsterbock 19891 ). How- 
ever, dust will absorb H/3 more efficiently than Ha, thus any dust present will 
increase this value. Comparing the measured value of Ha/H/3 to the theoretical 
value allows one to measure the amount of dust attenuation. However, Ha mea- 
surements are currently restricted to z < 2.7 for practical reasons, as at higher 
redshifts this line shifts out of the K-band and is out of reach from the ground. 
Currently, Ha measurements of high-redshift LAEs exist for only a handful of 
objects from the Double-Blind survey (Hayes et al. 2009). Thus, to truly probe 
the physical properties of LAEs, we need a large sample at z < 2.7. 
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Figure 3. The ratio of Lya/Ha versus the measured dust extinction from 
the Balmer decrement fo r 11 LAEs at z ~ 0.3, confirmed to be star-forming 
(jFinkelstein et al.ll2009al fbl). The different curves show how Lya/Ha changes 
with increasing dust extinction levels for various values of the ISM dumpiness 
parameter q. Thus, it appears that at least for these low-redshift LAEs, the 
inhomogeneity of the ISM may be enhancing their Lya EWs over the intrinsic 
value, as they are consistent with g < 1 (from Finkelstein et al. 2010, in prep). 
NIR spectroscopy will allow similar investigations into LAEs at z ~ 2-3. 



Beginning in the fal l of 2011 , the Hobby Eberly Telescope Dark Energy 



Experiment fHETDEX. iHill et al.l 1200811 wih provide these objects. HETDEX 
will discover nearly one million LAEs at 1.9 < z < 3.5 over ~ 500 deg^, in order 
to measure baryonic acoustic oscillations, and thus probe the evolution of the 
dark energy equation of state. Measuring the rest-frame optical emission lines 
from the LAEs at the lower end of this redshift range (possible in a few hours 
with 8-lOm class telescopes) will provide some of the first direct measurements 
of the dust extinction and the metallicities of LAEs. Additionally, via ratios 
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of rest-frame optical emission lines ( Baldwin. Phillips. Sz Terlevich 198ll ). one 
can also probe for AGN contamination, as low-luminosity AGN may be lurk- 
ing in the shadows. Finally, by comparing the ratio of Lya to a n on-resonant 
line ( i.e. Ha), one can di rectly probe the geometry of the ISM (cf. lAtek et al 



20091 : IScarlata et al.l[2009l '). Figure [3] shows 11 GALEX-discov ered LAEs at 



0.3 wh ich were determined to be star-formation dominated (jFinkelstein et al 



2009bl ) . In a uniform ISM, the ratio of Lya/Ha should decrease rapidly with 



increasing dust extinction, as shown by the q = 10 curve. However, as the 
ISM becomes more clumpy, Lya becomes less attenuated with respect to Ha 
(and the continuum), as shown by the curves for decreasing values of q, until at 
very small values of q, th e ratio of Lya /H a becomes greater than the intrinsic, 
theoretical value of ~ 9 ( Qsterboc^ Il989l ) . Using optical spectroscopy of these 
low-redshift LAEs, we determined that many of them appear to have ISMs which 
enhance their observed Lya EWs (Finkelstei n et al. 2010, in prep), consistent 
with their stellar population modeling results (jFinkelstein et all 2009al ) . At high 
redshift we do not currently have any direct probes into the state of the ISM, 
but rest-frame optical emission lines will provide these probes in the near future. 



4. Approaching the Big Bang 

The future of galaxy evolution studies diverges into two extremes. As described 
above, we need to study LAEs (and all galaxy types) at z < 3 in order to di- 
rectly measure their physical properties via diagnostic emission lines. At the 
same time, pushing galaxy studies to the highest redshifts possible can tell us 
more about how the first galaxies formed, and their impact on the last ma- 
jor phase transition of the Universe, the reionization of the IGM. Perhaps the 
best chance of discovering galaxies at z > 7 rests with the narrowband selec- 
tion technique, as the bright Lya emission can compensate for the very faint 
continuum levels. At these redshifts, Lya is observed in the NIR, thus strong 
atmospheric OH emission causes difficulties. The typical way to compensate for 
this is to manufacture ultra-narrowband filters (FWHM 10-100 A) which are 
sensitive to regions between these sky lines. The most promising gaps are at 
~ 1.06 and 1.19 /im, corresponding to Lya at z = 7.7 and 8.8. Many ongo- 



ing o r planned sur veys will exploit the se gaps, such as DAzLE (iHorton et al . 
20041 ) and ELVIS (jNilsson et al.1 l2007al ). Currentlv. iHibon et al.l (f2009l) found 



7 z ~ 7.7 LAE candidates down to fj^ya ~ 8 x 10~^^erg s~^ cm~^ with a 40 
hour 1.06 fim narrowband exposure with the CFHT. Using a 28 hour exposure 
obtained with the NEWFIRM wide-field NIR camera on the KPNO 4m, Tilvi et 
al. (2010, in p rep) discovered ano ther few candidate z ~ 7.7 LAEs. However, all 
but one (from iHibon et al.l l2009l ) of these candidates were only detected in the 



narrowband and not in any broadbands, although this is not unexpected given 
the limitations of ground-based NIR imaging. Thus, we require spectroscopic 
follow-up to determine the validity of these candidates. 

At z ~ 3-4, it is straightforward to obtain large samples of LBGs and 
LAEs. At z ~ 5-6, the dropping observed continuum brightness makes the LAE 
selection more powerful. However, at z > 7, we've gotten to the point where 
ground-based narrowband selection is becoming extremely difficult. Due to the 
recently installed Wide Field Camera 3 (WFC3) on board HST, we now have 
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Figure 4. 5" cu touts of a Zphot = 7.55 galaxy discovered with WFC3, from 
iFinkelstein et all (|2009d) - The right-hand plot shows the probability distri- 
bution function (PDF) of the photometric redshift. This object has > 99% of 
its PDF at z > 6. 



the ability to probe down to mAB ~ 29 in the NIR with broadband photome- 
try, thus LBG selection at the highest redshifts is now easier than narrowband 
LAE selection. Soon after WFC3 was commissioned, it obtained 60 orbits in 
the HUDF in th e Y, J and i J-baii ds, a t 1.05, 1.25 and 1.6 j^m , respectively 
(PI Illingworth). iQesch et alJ (l2009l ) and lBouwens et alJ (I2009bl 1 pubhshed the 
discoveries of 15 z ~ 7 (z' dropout) and 5 z ~ 8 (y dropout) L BG candidates. 



Amon g their conclusions were that these galaxies were very blue (jBouwens et al 



2009al ). Does this imply that we are finally observing the first galaxies? 



To take this one step further, we used an updated reduction of th e WFC3 



data to perform an independent photometric redshift a nalysis (similar to McLure et al. 
20091 ) ■ finding 35 candidate galaxies at 6.3 < z < 8.6 ( Finkelstein et al. 2009dl ). 



Figure m shows one of our galaxy candidates at Zphot = 7.55lQ'3g. We examined 
the rest-frame UV colors of t hese galaxi es, comparing them to the empirical lo- 
cal starburst templates of iKinnev et al.l (|l996 ) as well as a number of synthetic 
stellar populations using updated models from Bruzual & Chariot (2003). Fig- 
ure shows the J — H colors of the subsample of galaxies at z < 7.5 versus 
their ff-band magnitudes. Given that the photometric errors on the individual 
objects are rather large, we computed an average color for all objects (circle), 
objects with H < 28.5 (triangle) and objects with H > 28.5 (large square). 
These means and their corresponding errors were computed via bootstrap sim- 
ulations, where in each simulation we varied the color of the individual points 
by an amount proportional to their photometric scatter (which we determined 
with separate simulations), and then recomputed the mean color. 

It is immediately apparent that both the bright and the faint galaxies 
are significantly (> 4 (t) bluer than local starbursts. This is interesting, as 
Papovich et al.l ( 200ll ) performed a similar analysis on LBGs at z ~ 3, and 



found their colors to be consistent with the Kinney et al. templates. Rather, 
these galaxies are consistent with NGC1705, which is one of the bluest local star- 
burst known. The stellar population models, shown as the shaded arrows on the 
right-hand side, imply that these bluest galaxies are consistent with young ages 
(< 100 Myr) sub-solar metallicities and little-to-no dust extinction, representing 
a change from the more evolved LBGs typical at z < 6. 

To learn in more detail about the physical properties of these highest red- 
shift objects, we performed SED fitting using their observed WFC3 fluxes, as 
well as limits from ACS and IRAC. With only a few detected data points, defini- 
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Figure 5. A co l or-mag nitude diagram for galaxies at 6.3 < z < 7.5, from 
iFinkelstein et al.l ()2009dl ). The individual points are the small gray squares. 
We computed the average colors for all objects (circle), objects with H < 
28.5 (triangle) and H > 28.5 (large square). While LBGs at z 3 have 
colors consistent with local starbursts, all of these bins are significantly bluer. 
The faint bin is bluer than even NGC1705, which is one of the bluest local 
starbursts known. This implies that these galaxies, especially the faint ones, 
are consistent with having young, low metallicity and relatively unextincted 
stellar populations. 




Figure 6. The best-fit and maximum masses for our sample at z ~ 7 (left) 
and z 8 (right), showing a typical stellar mass of 10^-10^ Mq. These 
galaxies are lower mass than z 3 LBGs (shown by the arrow), and we even 
see a hint of evolution from z ~ 7, as z ~ 8 LBGs are more likely to have 
M ~ IO^Mr. 
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tive determinations of their ages, dust extinction and metallicities are difficult. 
However, they do appear to be young, consistent with the color analysis above, 
as 30/35 objects have very young ages of < lOMyr within the 68% confidence 
range on their age (as determined from 10^ bootstrap simulations). While the 
mean extinction is Ay ~ 0.4 mag, over half of the sample is consistent with no 
dust extinction. Lastly, while the metallicity is traditionally difficult to decipher 
in SED-fitting analysis, the very blue colors of these objects allow us to place a 
significant limit on the metallicity, with Z < O.IZq significant at the la level. 



1 1 F 




- OMqss of L*^alaxf 

2 4 6 8 

Redshift 

Figure 7. The stellar masses of L* galaxies versus redshift from Finkelstein 
et al. (2009d), with comparison samples taken from the literature. Our re- 
sults are the large stars at z ^ 7 and z ^ 8. The remainin g sy mbols, from 
dark to light, are: iReddv et all (|2006D . IShaplev et all (|2Q05[ ) and IStark et~all 
([2009) at z = 4, 5 and 6, respectively. The circles denote the mass at 
L* (charac teristic luminos i ty) at each redshift, using the lumi nosity func- 
tions from iBouwens et"aI1 (|2007[ ). iReddv fc Steidell (|2009t ) and lOesch et al.l 
(|20Q9[ ). The background gray circles de note stellar masses of Lyg emittin g 
galaxies at 3.1 < z < 6.5, taken from: Chary. Stern, & Eisenhardt ('2005'), 
iGawiser etall (I2006D. [Pi rzkal et al. (2007), Nilsson et al. (2007b), Lai et al] 
(|2007l ). lFinkelstein et all (|2007l ). lLai et al.l(|2008l ) fc lFinkelstehi'et'al.l (|2009e). 
The gray hatched region denotes the interquartile range of the LAE masses. 
The masses of the z > 6.3 LBGs studied here are more similar to those of 
LAEs at all redshifts than LBGs at any redshift < 6. 

While the above properties are difficult to directly pin down as they all result 
in reddening of the SED, constraints on the stellar masses are more robust, as 
the degenerate combination of all other parameters yields similar mass-to-light 
ratios. Figure [6] displays the distribution of best-fit masses at z ~ 7 and z ~ 8, as 
well as the cumulative probability distribution of stellar masses at each redshift. 
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This figure shows that the galaxies in our sample have typical masses of 10^- 
10^ Mq, with galaxies at 2 ~ 8 having a higher probability of having masses 
of < 10^ Mq. However, these masses were derived using solely rest-frame UV 
fluxes, which can miss significant amounts of mass in older stars lost beneath the 
glare of the UV-bright young stars. We thus performed a second fit, allowing two 
bursts of star formation, where one burst was constrained to occur at z = 20, 
and contain 90% of the stellar mass. We refer to the best-fit masses from these 
fits as the maximum mass, as they contain the maximum physically plausible 
amount of mass which could exist in these galaxies without violating the Spitzer 
flux limits (although physically unlikely, we did perform a fit with 99% of the 
mass in old stars, and found the maximum mass to only increase by a factor of 
~ 2). Investigating Figure [6l one can see that these maximum masses are only 
shifted to marginally higher masses than the best-fit masses, implying that the 
Spitzer upper limits are providing little constraints on the stellar masses. This 
is understandable, as at these redshifts, the time elapsed since 2; = 20 is only 
~ 400-700 Myr, thus late- type B/early-type A stars will still be on the main 
sequence, and these stars contribute measurably to the rest-frame UV. 

Figure \7\ shows the masses of L * galaxies as a furi c tion of redshift, with 
the lower redshif t poiri ts coming from lReddv et al.1 (I2006l 1. IShaolev et al.1 (l2005l ) 
and IStark et al.l ( 20091 ). Investigating this figure, from z = 2~6 the mass of a 
typical L* galaxy appears to drop by an order of magnitude, though the spread 
in the z ~ 6 point made an exact determination difficult. With the addition of 
our data points, we now confirm this trend of decreasing mass with increasing 
redshift, with a slight (though not significant) drop in mass from z ~ 7-8. The 
small gray circles denote the masses of LAEs taken from the literature, and the 
hatched box denotes their interquartile range. At 2; < 5, LBGs are systematically 
more massive than LAEs. However, at 2; > 6, there is a change in the physical 
properties of typical continuum selected galaxies, and their masses are more 
consistent with LAEs at all redshifts, than with LBGs at any previous redshift. 
Given the blue colors and corresponding un-evolved characteristics of the z ~ 7- 
8 galaxies discussed here, it appears that at these high redshifts, evolved LBGs 
common at lower redshifts no longer exist, and our sample consists of galaxies 
which are the true building blocks of lower-redshift evolved galaxies, similar to 
LAEs at lower redshifts (e.g.. .Gawiser et al..,2007i ). 



5. Conclusions 



Selecting galaxies on the basis of their Lya emission allows the discovery of a 
plethora of faint, high-redshift galaxies. Numerous studies over the past few 
years have shown that these galaxies primarily appear to be a unique popu- 
lation of star-forming galaxies, which are young (< 100 Myr) and low mass 
(< W^Mq), although a minority can be older and more massive. Many of the 
LAEs which have been discovered exhibit strong Lya equivalent widths which 
cannot be explained by normal stellar populations. While AGN contamination 
could reconcile this problem, AGN (at least luminous ones) appear to be rare 
among LAEs. Very low metallicities could also cause high EWs, but the dusti- 
ness of many LAEs likely rules out this scenario. Finally, this dust, if in a 
clumpy geometry, can enhance the Lya EW by increasing the escape fraction 
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of Lya relative to the continuum, and this scenario can explain the SEDs of 
a sample of high-redshift LAEs. However, to truly determine the make-up of 
LAEs, including their dust properties, we require direct measurements, which is 
available presently with NIR spectroscopy, and in the near future with ALMA. 

Pushing galaxy studies to the highest redshifts, ground-based narrowband 
selection of LAEs becomes extremely difficult due to the bright NIR sky back- 
ground, with the deepest studies finding only a handful of candidate z > 7 LAEs. 
However, the addition of WFC3 to HST provides the capability to obtain large 
samples of z > 7 galaxy candidates. Initial studies of these galaxies have shown 
them to be somewhat primitive in their make-up, with young ages (< 100 Myr), 
low masses (10^-10^ ^-^0)1 low mctallicity (< O.IZq), and consistent with little- 
to-no dust extinction. Comparing the physical properties of these most distant 
galaxies to those at lower redshifts, they appear more similar to LAEs at all 
redshifts than LBGs at any redshift. This implies that typical galaxies at z > 7 
are less evolved than at lower redshifts. At some point in the past we should 
begin to see the first building blocks of galaxies, and these observations show 
that z =7-8 may be this epoch, although larger samples of spectroscopically 
confirmed galaxies will be necessary to make a definitive determination. 
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